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Abstract. The age-metallicity relation of the solar 
neighbourhood is studied using a sample of 552 late-type 
dwarfs. This sample was built from the intersection of pho- 
tometric catalogues with chromospheric activity surveys 
of the Mount Wilson group. For these stars, metallici- 
ties were estimated from uvby data, and ages were calcu- 
lated from their chromospheric emission levels using a new 
metallicity-dependent chromospheric activity-age relation 
developed by Rocha-Pinto & Maciel ( fl998 ). A careful es- 
timate of the errors in the chromospheric age is made. The 
errors in the chromospheric indices are shown to include 
partially the effects of the stellar magnetic cycles, although 
a detailed treatment of this error is still beyond our knowl- 
edge. It is shown that the results are not affected by the 
presence of unresolved binaries in the sample. We derive 
an age-metallicity relation which confirms the mean trend 
found by previous workers. The mean metallicity shows a 
slow, steady increase with time, amounting at least 0.56 
dex in 15 Gyr. The initial metallicity of the disk is around 
—0.70 dex, in agreement with the G dw arf metallicity dis- 
tribution (Rocha-Pinto & Maciel 1996 ). According to our 
data, the intrinsic cosmic dispersion in metal abundances 
is around 0.13 dex. a factor of two smaller than that found 



by Edvardsson et al. ( 1992 ). We show that chromospheric 
ages are compatible with isochrone ages, within the ex- 
pected errors, so that the difference in the scatter cannot 
be caused by the accuracy of our ages and metallicities. 
This reinforces some suggestions that the Edvarsson et 
al.'s sample is not suitable to the determination of the 
age-metallicity relation. 

Key words: stars: late-type - stars: chromospheres - 
Galaxy: evolution - solar neighbourhood 



1. Introduction 

In the general picture for the evolution of our Galaxy, the 
stars form from gas enriched by previous stellar gener- 
ations, and eject new metals to the interstellar medium 
after their death. Due to this continuing enrichment, old 
stars are likely to be metal-poorer than younger stars. 

The age-metallicity relation (hereafter AMR) has been 
the subject of many studies in the literature. The first 
systematic attem pt to determine this relation was made 
by Twarog (198C). He presented photometric metallicities 
and isochrone ages for 329 F dwarfs, finding a smooth 
increasing relation with an average scatter of 0.12 dex. 
Twarog's sample was reanaly sed su bsequently by two dif- 
ferent groups. Carlberg et al. (|1985|) have found a very flat 
AMR, probably because they have cut from the sample all 
stars with metallicities lower than —0.50 dex. On the other 



hand, Mcusinger et al. ( 1991 ) used updated isochrones and 
a metallicity calibration and found an AMR very similar 
to that of Twarog. Other attempts to derive this relation 
using photometric metallicities were made by Ann & Kang 
( 1985 ) and Marsakov et al. ( |199C| ) , but these works suffer 
from the lack of an unbiased sample selection. 

By far the most common approach to study the AMR 
is the use of photometric metallicities and isochrone ages, 
since this allows for a large sample which can compensate 
for a poor accuracy in these quantities. However, some 
stu dies m ake use of spectroscopic metallicities. Nissen et 
al. ( 1985 ) have presented [Fe/H] and ages for 29 F dwarfs 
taken from the larger sample that was investigated later in 
more detail by Edvardsson et al. ( 1993 , hereafter Edv93). 
Their data agrees well with Twarog's AMR, although the 
scatter seems to be higher. Lee et al. (1989) give a spec- 
troscopic AMR for 559 disk stars whose metallicities are 
given in the catalogue of Cayrel de Strobel et al. ( 1985 ). 
The authors measure ages from isochrones in five different 
diagrams, some of which are likely to be independent of 
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error in the isochrone age determination (see for exam- 
ple Ng & Bertelli 1998| ). However, their sample is highly 
heterogeneous, not only regarding the metallicity sources, 
but also the spectral types used for the study. 

Presently, the most significant work on the AMR was 
done by Edv93. They measured accurate spectroscopic 
metallicities on 189 carefully selected disk stars. Ages were 
found by VandenBerg's (1985) isochrones. Their result 



is rather surprising: while the mean AMR is similar to 
Twarog's AMR, the metallicity dispersion is so high that it 
casts doubts about the real meaning of the age metallicity 
relation. The s ame data were recently reanalysed by Ng & 
Bertelli ( 1998 ), using updated isochrones as well as HIP- 
PARCOS parallaxes, and the high dispersion in metallicity 
was confirmed. A high metallicity dispersion can also be 
found i n the AMR of open clu sters (Strobel 1991 ; Jane s 
& Friel |1993| ; Carraro & Chiosi |1994| ; Carraro et al. |1998| ). 



Although this seemed puzzling with respect to the pre- 
vious well-marked photometric AMR, it has given rise to 
a new era in studies about the chemical evolution of the 
disk in which the gas is not very well mixed so that local 
inhomogeneities would obscure the overall growth of the 
metallicity. 

Nevertheless, the AMR is still a poorly known func- 
tion. A critical review of the literature shows that only two 
independent carefully selected samples of field disk stars 
were ever used in its study: Twarog's and Edv93's sam- 
ples. And the conclusions that can be drawn from both 
samples, mainly on the metallicity dispersion, are in re- 
markable disagreement, strongly motivating the present 
study. One other aspect that deserves further investiga- 
tion is the apparent lack of metal rich objects in t he ag e 
range 3-5 Gyr as pointed out by Carraro et al. ( 1998[ ), 
since a similar feature is als o marginally visible in the 
data of Carlberg et al. ( 1985 ). 

This work is the first of a series in which we revisit 
some fundamental constraints on chemical evolution, by 
using an independent and extensive sample of long-lived 
dwarfs. The main novelty of these papers is an attempt 
to explore an independent tool to measure stellar ages: 
the chromospheric activity level, which can be even more 
accurate to measure ages of late-type stars than isochrones 
(Lachaume et al. 1999 ). The only study ever published 
using this approach is that of Barry (1988) but the sample 
selection prevents his AMR to be taken as representative. 
There are several chromospheric indices in the literature. 
Here, we will use both S and logi?Jj K indices, as defined 
by Noyes et al. (1984). The first is a measure of the line 
intensity in the H and K Ca II related to the continuum 
(Baliunas et al. 1995). From its definition, S depends upon 
the stellar colour. Noyes et al. provide equations for the 



HK' 



which is a colour- 



transformation of this into 
independent index. 

The paper is organized as follows: section 2 describes in 
detail the data, where we pay special attention to the esti- 



sample. In section 3, the AMR is derived and a compari- 
son with previous relations in the literature is presented. 
A detailed discussion about the magnitude of the metal- 
licity dispersion is given in section 4. The final conclusions 
follow in section 5. 

2. The sample 

2.1. Selection criteria 

The criteria we have followed for the construction of the 
sample were based on the requirement to have a number, 
as large as possible, of disk stars for which reliable metal- 
licities and ages could be determined. This was done by 
taking the common stars in the photometric catalogues 
of Olsen ( |1983| , |l993j , |1994|) and the chromospheric activ- 
ity surveys of Soderblom (1985, hereafter S85) and Henry 
et al. ( |1996[ hereafter HSDB). This procedure yielded 
an initial sample composed by 729 late-type dwarfs with 
0.307 < (b-y) < 0.622 and -5.40 < \ogR' KK < -3.78. 
Here, (b — y), as well a s mi and c\ refer to the standard 
uvby indices (Crawford 1975 ). 

According to HSDB, the late- type dwarfs in the solar 
neighbourhood can be divided according to four chromo- 
spheric populations, namely the very active stars (logi?Jj K 
> -4.20), the active stars (-4.20 > logi?^ K > -4.75), 
the inactive stars (—4.75 > logi?^ K > —5.10) and the 
very inactive stars (logR' HK < —5.10). We will refer to 
these groups, hereafter, as VAS, AS, IS and VIS, respec- 
tively. 

In Figure [l]a,, the diagram (b — y) x c\ for these 729 
stars is shown. The polygon in this figure corresponds to 
the area of the dia gram occupied mostly by subgiants, 
according to Olsen (1984). As much as 11% of the sam- 
ple is probably composed of subgiants. Figure [l]b shows a 
histogram of chromospheric activity levels for these stars. 
Many of them have logR^ < —5.00, corresponding to a 
chromospheric age greater than 5.6 Gyr, according to the 
age calibration by Soderblom et al. (1991), an additional 
evidence for their evolved status. 

One of these supposed subgiants (HD 119022, shown 
in Figure |l|a as an open triangle) is a VAS (logi?^ K = 
—4.03), and presumably should be a very young star, not 
a subgiant. Note that the star is quite near the limit of the 
subgiant area, and could rather be considered a slightly 
evolved dwarf. Its nature is still uncertain, since it is lo- 
cated far be yond the zero-age main sequence. Soderblom 
et al. ( 1998 ) have studied this star in more detail, and have 
found spectral features typical of youth (strong Li and 
broad- line spectrum), but its high luminosity prevented 
them from classifying it unambiguously as a young star. 
Other trends in its spectrum (strong Ha in absorption, a 
discrepancy between colour and spectral type, and deep 
sharp features inside the Na D profiles) suggest that it 
could be a heavily reddened star. We will return to this 
star later, and show that the reddening does affect it sub- 
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Fig. 1. Photometric and chromospheric data for the ini- 
tial sample of 729 late-type dwarfs, a) Diagram (b — y)x c\ 
for the initial sample. The polygon indicates the area pop- 
ulated by subgiants, and the open triangle shows the po- 
sition of the star HD 119022 in this diagram (see text), b) 
Distribution of chromospheric activity indices of the stars 
located in the subgiants area. 



These 'subgiant' stars are probably a mix of thick disk 
and old thin disk stars. We choose to keep all of them in 
the initial sample since their number is relatively small and 
because they most probably represent the older epochs of 
the disk evolution we are trying to recover. However, it 
is not presently known whether their chromospheric prop- 
erties can be related to that of main-sequence stars. At 
least, as will be demonstrated in the next pages, there is 
good agreement between their isochrone ages and chromo- 
spheric ages. 



2.2. Metallicities and cole 



idices 



The calibrations of Schuster & Nissen ( |1989| ) were used 
for the determination of the metallicity of each star. The 
metallicities for 3 stars redder th an (b — y) = 0.599 were 



To obtain the colours 5m\ and Sci, we have adopted the 
standard c urves (b — y) x m\ and (b — y) x c\ given by 
Crawf ord ( 1975 ) for late F and GO dwarfs, and by Olsen 
( 1984 ) for mid and late G dwarfs. The error in the metal- 
licity from these calibrations is expected to be around 0.16 
dex (Schuster & Nissen |1989|) . 

To account for the mi deficiency in active stars (Gime- 
nez et al. 1991, and references therein), the photometric 



metallicity of the active and very active stars (log i?^ K > 
—4.75) needs to be corrected by adding to it an amount 
A[Fe/H], given by the equation proposed by Rocha-Pinto 
& Maciel (|1998[): 



A[Fe/H] = 2.613 + 0.550 log R' HK . 



(1) 



Since this correction is still rather uncertain, we will 
always present the uncorrected metallicity (which is the 
quantity that comes directly from the photometric data) 
in all plots, tables and references in the text, except when 
mentioned otherwise. 

2.3. Chromospheric ages 

Chromospheric ages were calculated using the age calibra- 
tion developed by Soderblom et al. ( 1991 ; their Equation 
3). Besides this calibration, we have tested the calibration 
by Donahue (1993, as quoted by HSDB). This last calibra- 
tion gives ages very similar to those of Soderblom et al. 
( 1991 ), except for stars with logi?Jj K > —4.40, for which 
it predicts ages systematically lower. We have not used 
it in our final results for the sake of consistency, as some 
corrections applied to the chromospheric ages were based 
on the calibration by Soderblom et al. ( |l99l| ). Moreover, 
only 6% of the sample stars have log R[ 



HK 



> 



-4.40. We 



have also verified that our conclusions are not dependent 
on the difference between these age calibrations. 

The adopted age calibration was further corrected for 
the dependence on the metallicity proposed by Rocha- 
Pinto & Maciel (199S). In Figure we show a compari- 
son between the uncorrected age calibration, and the same 
curve corrected for three values of [Fe/H], namely —1.2, 
—0.23 and +0.4, which corresponds to the lowest, the 
average and the highest metallicity of our sample. Also 
presented in the same plot are the age calibrations by 
Soderblom et al. (1991; their Equation 1) and that by 
Donahue. Most important in this plot is that it shows 
that the majority of our stars will have an age not very 
differe nt fro m that given by the uncorrected Soderblom 
et al. (1991)'s age calibration. This can be seen from the 
close agreement between the uncorrected age calibration 
and that for the average metallicity of the sample. 

Figure || shows the resulting age distribution for all 
these stars, excluding 54 stars that present ages greater 
than 15 Gyr. The unrealistically high age of these 54 stars 
could be caused by one of these reasons: (i) the star is 
experiencing a Maunder-minimum phase (see Baliunas et 
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Fig. 2. Comparison between chromospheric age calibra- 
tions in the literature. Solid line, Soderblom et al. (1991 



equation 3); circles, Soderblom et al. (1991, equation 1); 
stars, Donahue (1993); dashed line, adopted calibration 
for three metallicities, namely —1.2, —0.23 and +0.4. 




Fig. 3. Age distribution for the initial sample. Stars with 
chromospheric ages greater than 15 Gyr are not shown in 
this plot. 



[Fe/H] worked together to produce a spuriosly high age; 
and (iii) the age calibration could be overestimated for 
large log R^ K values, since at that range it is given by an 
extrapolation of data for younger stars. 

Several errors contribute to the error in age, coming 
from the procedures of corrections and transformation of 
indices into age. We consider two error sources, namely the 
error in the index log i?HK an d the error in the photomet- 
ric metallicity which enters in the mctallicity-dependent 
correction. Rigorously speaking we should also consider 
the error in Soderblom et al. (1991)'s calibration and the 
error in the metallicity correction itself. However, these 
are not independent from the former, since the scatter in 
the calibration reflects mainly the neglecting of the metal- 



The average error in log -Rh k was estimated from the 
data published by Duncan et al. ( 199l| ). They present data 
for the variation of the S index in several late-type stars, 
for a time interval of 17 years. We have calculated the 
corresponding values for log i?HK using the equations by 
Noyes et al. ( [1984 ). The index shows much variation for 
some stars, but the average error is around ±0.05. This er- 
ror estimate agrees closely with that made by S85. The un- 
certainty introduced in the age by the photometric metal- 
licity was estimated by propagating the error in [Fe/H], 
using the metallicity-corrcctcd age calibration given by 
Rocha-Pinto & Maciel (1998). 

A further complication arises if we consider other error 
sources: flaring, rotational modulation of active regions, 
and the stellar magnetic cycles. A small part of this error 
is already incorporated into the error in log R'n K , but some 
stars show much variation in the chromospheric indices. 
Donahue ( 1998|) shows that the age of the Sun could be 
miscalculated within an error of several Gyr if its chromo- 
spheric age were derived from observations in an epoch of 
maximum or minimum activity. S85 remarks that the most 
serious problem is with the stellar magnetic cycles, and the 
effects of the other sources of variability are smaller. This 
uncertainty is the major problem of the chromospheric 
ages. It can be much decreased provided that we have a 
good determination of the mean stellar activity, which is 
still not the case for the majority of the stars in the solar 
neighbourhood. The promising accuracy of this technique 
was presented also by Donahue (1998), who has shown 
that the age discrepancy between binaries is usually lower 
than 0.5 Gyr for systems younger than 2 Gyr, and does 
not exceed 1.0 Gyr for older pairs. 

We have investigated this error with more detail in the 
data published by Duncan et al. ( 199l| ), for 85 stars that 
have been continuously monitored from 1966 to 1983. The 
logR' HK index for each of these stars was calculated, and 
an average error was found. The results are presented in 
Figure |4|, where we have separated the stars according to 
their colours. The Sun is also shown in this diagram (in 
the bottom-left panel). To calculate its position, we have 
used the S indices for the maximum and minimum activity 



during cycles 20-22, given by Donahue (1998). 

Note that there is a small tendency to find the biggest 
errors in the redder stars. This pattern follows more or less 
closely the findings by Baliunas et al. (1995), who have 
shown that F dwarfs do not show much magnetic vari- 
ability, while the redder dwarfs present very well defined 
magnetic cycles. This can increase the error in log R'n K 
beyond the value we have adopted, since from only one 
observation it is impossible to know exactly in what part 
of the cycle the star is. On the other hand, it is evident 
from the figure that no single law can be used to estimate 
the error in the index, given the chromospheric activity 
level and the stellar colour. 

The \ogR' HK error for the Sun is one of the greatest. It 
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Fig. 4. Error in t he ind ex logi?[j K due to the stellar magnetic cycles, calculated from the data for 85 stars published 
by Duncan et al. ( 1991 ). A small trend is seen in the sense to find the greatest errors in the old redder dwarfs, but the 
behaviour is not so well defined to allow the use of a formalism describing it. The Sun is also shown in this diagram, 
in the bottom-left panel. 



knowledge of the solar magnetic activity (for instance, the 
time span for the solar S measurements in this figure is 
about 3 times longer than that of the stars), or to a real 
more pronounced magnetic variability in the Sun. 

Around 80% of our sample is composed of stars having 
only one log -R^k measure. Their ages are most subject to 
the errors due to stellar magnetic cycles. Nevertheless, the 
incorporation of this kind of error is very difficult, and we 
have decided to use a conservative value of 0.05 dex cal- 
culated as described above. This error is 0.01 dex greater 
than that estimated by S85 for 33 dwarfs. According to 
this author, on the average the star will present an un- 
certainty around 10% in \ogR^ K (0.04 dex), due to the 
stellar magnetic cycles. This procedure will not affect sig- 
nificantly our conclusions, since the error in log R'n K is 
small compared to that due to the photometric metallic- 
ity. 

The impact of the individual error sources on the age is 
shown in Figure ||[ which shows relative errors (e ago /age), 



line, was calculated by adding those two individual errors 
in quadrature. This is the error estimate used throughout 
this paper. As can be seen from the error trends in the 
plot, the positive error in the age is greater than the neg- 
ative error. This is due to the logarithmic nature of the age 
calibration and the metallicity corrections. We expect that 
many stars will show ages scattered over a large range of 
values above the real stellar age. This is one of the reasons 
why some stars present unreasonable ages greater than 15 
Gyr. 

Since the chromospheric age depends not only on the 
activity level of the star but also on its metallicity, the 
task of finding a chromospheric age is analogous to find 
an 'isochrone' age in the metallicity-activity diagram, as 
shown in Figure |^. This figure shows lines of equal chromo- 
spheric age, expressed in Gyr. The dotted vertical lines are 
used to separate the four chromospheric populations, as 
defined section 2.1. An interesting result is that while for G 
dwarfs the younger isochrones are crowded in the colour- 
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in the index log B! UK (dotted lines), and in the photometric 
metallicity (dashed lines). The final age error was calcu- 
lated by the root mean square of all these errors and is 
shown as thick solid lines. 




Fig. 6. Curves of equal chromospheric age in the 
metallicity-activity diagram. The average error in this dia- 
gram is shown at the lower left corner. The vertical dotted 
lines separate the four groups of stellar activity, as defined 
in section 2. 

rate isochrone ages for very young stars, they are much 
more spaced in the metallicity-activity diagram. The op- 
posite occurs for the older isochrones being more crowded 
in this last diagram. From this we can see that chromo- 
spheric ages are most useful for younger stars, losing part 
of their accuracy as we consider older stars. 

2-4- Parallaxes, absolute magnitudes and reddening 

Parallaxes and absolute magnitudes were obtained from 
the HIPPARCOS data for 714 of these stars. The his- 
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Fig. 7. Histogram of distances from the Sun, for the initial 
sample. The majority of the stars are located within 60 pc, 
but 5% of them have distances greater than 100 pc. 



the stars from the Sun. The majority of the stars are lo- 
cated within 60 pc, but 5% of them have distances greater 
than 100 pc. 

We have calculated the reddening for 13 of these stars 
located bey ond 1 00 pc, using (3 indices from Hauck & 
Mermilliod ( 199E ) a nd th e intrinsic colour calibration of 
Schuster & Nissen ( 1989 ). The results are presented in 
Table [l] where we list the HD number, the distance from 
the Sun, log-R^K: the colour excess E(b — y), the dered- 
dened Stromgren indices (b — y)o, mo and Cq, the photo- 
metric metallicity calculated with the dereddened indices, 
and the difference between the previous [Fe/H] calculated 
without and with reddening correction. Four of these stars 
present E(b — y) < and should not be reddened. Thus, 
dereddened indices are not provided for them in the Table. 

The Table shows that there is considerable reddening 
for some stars, and this can seriously affect the stellar 
metallicity. Provided that we had (3 for all the stars lo- 
cated beyond 100 pc, we could deredden their indices and 
find more reliable estimates for [Fe/H], and consequently 
age, but this index is available only for a few stars. To 
avoid reddened stars, in the next sections we will consider 
generally only those stars located within 80 p c, for which 
reddening is expected to be negligible (Olsen 1984 ). 

Note that the reddening for HD 119022 is apprecia- 
ble, confirming the suggestion by Soderblom et al. (1998). 
Moreover, taking the dereddened colours in Table |1|, we 
see that the star resides outside the subgiant polygon in 
Figure [l]a, right in the bulk of main-sequence stars. This 
can be taken as additional evidence for a lower age. 

Five objects amongst the 13 VAS of our sample have 
distances greater than 100 pc. We expect that they are 
mildly to strongly reddened. The straightforward con- 
sequence of this is that their photometric metallicities 
would seem lower than they really are. This can be 
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Name 


d(pc) 


log ^HK 


h{b - y) 


[o - y)o 


mo 


C() 


rriV. /TTl 

F e / H Jdcrcd 


Z\ |r e / rij 


HD 3611 


104.6 


— 5. 150 


0.014 


0.374 


0.166 


0.404 


-0.386 


-0.030 


HD 17169 


107.9 


— 4.720 


0.040 


0.432 


0.231 


0.332 


-0.044 


-0.300 


HD 39917 


194.9 


-4.050 


0.064 


0.454 


0.254 


0.298 


-0.180 


-0.444 


HD 119022 


124.2 


-4.030 


0.054 


0.420 


0.231 


0.336 


+0.053 


-0.397 


HD 122683 


176.7 


-4.760 


-0.018 












HD 139503 


109.9 


-5.180 


0.017 


0.376 


0.182 


0.385 


-0.057 


-0.144 


HD 141885 


104.1 


-5.340 


0.019 


0.390 


0.208 


0.459 


+0.155 


-0.160 


HD 142137 


105.8 


-5.040 


0.024 


0.381 


0.194 


0.435 


+0.059 


-0.210 


HD 151928 


114.8 


-5.070 


-0.023 












HD 159784 
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-5.200 


-0.032 












HD 179814 


133.3 


-4.970 


-0.017 












HD 199017 


173.6 


-4.920 


0.026 


0.412 


0.213 


0.325 


-0.072 


-0.195 


HD 202707 


112.0 


-5.110 


0.027 


0.396 


0.185 


0.343 


-0.188 


-0.225 



39917 and HD 119022, with a mctallicity difference of 
around —0.4 dex. Rocha-Pinto & Maciel ( 1998| ) showed 
that the photometric metallicity distribution of the VAS 
is peculiarly concentrated at lower [Fe/H], in strong con- 
trast with the expected youth of such stars. We have ex- 
plained this trend as the result of the mi deficiency, which 
makes stars with strong chromospheres resemble metal- 
poor stars from a photometric point of view. This defi- 
ciency can originate from the filling in of the line cores due 
to a chromosphere-driven photospheric activity (Basri, 
Wilcots k Stout |1989[ ). 

However, these new results suggest that the low pho- 
tometric metallicitics of the VAS could be explained as 
an effect of the reddening. At least, the difference in the 
metallicity, A [Fe/H], is large enough to make the photo- 
metric metallicity of the VAS very similar to the metallic- 
ity of the AS, in the [Fe/H] x logi?g K diagram. We need 
to know if the reddening can explain the 'activity strip' 
found in this diagram (Rocha-Pinto & Maciel 1998Q , since 
our procedure to correct [Fe/H] in the AS and VAS, for 
the effects of the ?»i deficiency, depends on the meaning 
of this strip. 

Only one other VAS has a published (3 in the litera- 
ture. It is HD 123732 which shows a very small reddening, 
E(b — y) = 0.002, in agreement with its distance from the 
Sun (around 63 pc). However, its photometric metallic- 
ity is not very low ([Fe/H] = —0.220) compared to other 
VAS. If we were to trust its youth, we would need to ex- 
plain its subsolar [Fe/H] as an evidence towa rds th e mi 
deficiency. Notwithstanding, Soderblom et al. ( 1998| ) have 
found features in its spectrum that classify it as a W UMa 
star, so that its chromospheric activity must not come 
from youth. On the other hand, Soderblom et al. measured 
spectroscopic [Fe/H] in two VAS which are most probably 
very young single stars, and the difference between these 
metallicities and the photometric metallicities are in very 



are located within 50 pc, and therefore are unreddened. 
The difference between the spectroscopic and photometric 
[Fe/H] can only be understood as resulting from the m\ 
deficiency. 

It is important to bear in mind that what Rocha-Pinto 
& Maciel ( 1998 ) have found seems to be a systematic trend 
depending on \ogR'^ K : the more active the star is, the 
more metal-deficient it looks. This trend makes the m\ de- 
ficiency hypothesis very appealing, since it can produce a 
similar effect. If this behaviour is likely to be produced in- 
stead by reddening, then we need to consider a correlation 
between distance and logf?y K . Figure || shows that such 
a correlation probably does not exist. The most distant 
VAS are not the most active. In fact, even disregarding 
the most distant VAS, the [Fe/H] x logi?[j K diagram still 
presents the 'activity strip'. The average metallicity of the 
VAS within and beyond 80 pc (our chosen distance cutoff) 
is —0.318 and —0.480 dex, respectively, showing that while 
reddening can explain part of the low photometric metal- 
content of the VAS, there is still another effect to account 
for, and the mi deficiency seems the most promising one. 

Since the proposed correction for the mi deficiency was 
empirically determined using only AS stars, and they do 
not seem to be affected by reddening, we have no reason 
to disregard them. Indeed, the agreement of our proposed 
correction with the two VAS studied by Soderblom et al. 



(1998) is good evidence that the extrapolation of the re- 
lation for the VAS is reasonable. 



2. 5. Effects of unresolved binarity 



According to Duquennoy & Mayor ( |l99l| ), 65% of the 
G dwarfs in the solar neighbourhood are presently un- 
resolved binaries. They are expected to present colours 
contaminated by the unresolved secondary, which would 
introduce errors in their metallicities and ages. We are 
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Fig. 8. Activity level versus distance from the Sun, for 
the VAS. The plot shows that a correlation between these 
quantities probably does not exist. The most distant VAS 
are not the most active. While reddening can explain part 
of the low photometric metal-content of some VAS, there 
is still another effect to account for, and the mi deficiency 
seems the most promising one. 
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Fig. 9. Expected number of binaries in each age bin, 
shown by dark bars. The dashed line shows the total num- 
ber of stars in each bin. The solid line indicates the ex- 
pected relative age error of the unresolved stars, as a func- 
tion of age. The error is always positive, since unresolved 
stars appear younger than they really are. 



more easily estimated. Since there is no set of chromo- 
spheric measurements for combined and resolved binaries, 
this discussion is predominantly theoretical. 

In principle, for two stars with the same age and metal- 
licity, logi?^ K are equal. The index that depends on the 
colour is the S index. For two coeval stars, born from 
the same parental cloud, S will be greater in the cooler 
one. The index is very dependent on the stellar colour, 
since the chromospheric flux is much greater in the late 



S index of the primary will be contaminated by that of 
the secondary. The combined index will be always greater 
than that of the primary, and the star would appear more 
active and younger than it really is. 

We have calculated the S index that would correspond 
to a certain logi?Jj K for several binary pairs, with (B — 
V) varying from 0.4 to 1.0, by inverting the equations 
provided by Noyes et al. (1984). A larger colour range 



could not be used due to the limit of the S calibration to 
colours higher than (B — V) = 1.0. The (B — V) colour was 
also used to estimate roughly the effective temperature 
and mass of the star. We thus can calculate the mass ratio 
for each pair (secondary to primary). 

To combine the S index of the two stars in the pair, 
we used a weighted mean. The weight is the Planck func- 
tion integrated from AA 3880 to 4020 A, which correspond 
to the spectral range where the Ca II lines are found. In 
this way, we can found an average S for the unresolved 
pair, and the corresponding \ogR' HK . We compared the 
difference between this average \ogR' HK and the real in- 
dex, known a priori. Basically, the difference is greater 
for the active stars, since the S index of the secondary 
becomes much greater than that of the primary. For the 
inactive stars, it is lower than 0.04 dex which is smaller 
than the error for the index presented in Figure ||. The 
difference depends also on the mass ratio, and on the pri- 
mary mass. A larger mass ratio tends to increase the dif- 
ference. That means that a pair of stars composed of a 
G+K dwarf will appear more active than that composed 
by two G dwarfs. But at a certain mass ratio, around 0.55, 
the behaviour changes and the difference diminishes. That 
is, a pair composed of a G+M dwarf would be much less 
contaminated than that composed of a G+K dwarf. This 
reflects the fact that the intensity of the spectrum of the 
hotter star becomes much more important than that of 
the fainter. The dependence on the primary mass is sim- 
pler: the more massive is the primary, the weaker is the 
contamination by the secondary, and the smaller is the 
difference in log -Rhk ■ 

Since these are just rough estimates, we decided to 
keep the qualitative approach, instead of deriving numeri- 
cal corrections. We have calculated the number of stars in 
each age bin that can have a wrong age due to unresolved 
binaries. This number is estimated as 



N b = ANf hi nf m<1 .Jq, 



(2) 



where AN is the number of stars in each age bin; /bin is 
the fraction of unresolved binaries in the sample, taken 
as 0.65; f m<1 _ 2 1S the fraction of stars with masses lower 
than 1.2M®, which are expected to be much influenced to 
contamination by a secondary, according to our calcula- 
tions; and f q is the percentage of mass ratios in which the 
primary S colour are affected by that of the secondary. 
From a Salpeter IMF, we estimate that f m<12 ~ 0.67. 
The fraction f q is conservatively estimated as 0.50 from 
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two stars with mass ratio around 1 will not affect the re- 
sults: they would have the same index; and stars much 
fainter will not contaminate the index of the primary. 

Figure [9] shows the number of stars with probably 
wrong ages due to duplicity. The solid line in this plot 
gives the expected relative age error for those stars. It uses 
the maximum error in logi?^ K , after varying the primary 
mass and the mass ratio. The error is always positive, since 
unresolved stars appear younger than they really are. 

It can be seen that the error is greater for the youngest 
stars. However, its magnitude is negligible in view of the 
age errors already considered in Figure |[ The number of 
stars subject to this errors is also small, so that we can 
conclude that it does not affect the main results of this 
paper. 

2.6. Representativeness of the sample 

An additional sample, consisting of 267 stars, was built 
in order to complement our initial sample. The sample 
consists of stars initially in the surveys of S85 and HSDB, 
but not having photometric data in Olsen's catalogues, as 
well as other stars scattered amongst several papers by 
the Mount Wilson group (Soderblom et al. 1991; Duncan 
et al. |1991| ; Baliunas et al. |1995| ; Saar & Donahue |1997| ), 
including the Sun itself which has not entered in the initial 
sample. In some cases, only the (S) index was provided, 
and we calculated the corresponding log R'n K indices using 
the formalism of Noyes et al. ( fl984 ). Photometric data 
for these stars were taken from the catalogue of Hauck & 
Mermilliod p9§). 

In figure |10| we show the main characteristics of both 
the initial and the additional sample: the metallicity and 
log i?HK distributions and the metallicity-activity dia- 
gram. In all three aspects considered, both samples differ 
somewhat. The metallicity distribution of the additional 
sample is broader, and the metallicity-activity diagram 
seems more scattered than that for the initial sample. Part 
of this scatter probably reflects the heterogeneity of the 
photometric data in the catalogue of Hauck & Mermilliod 
( |1998| ). Moreover, there are also significant differences be- 
tween the chromospheric activity distribution of the two 
samples. One of these differences is the excess of stars with 
—4.75 < logR' HK < —4.60, where the Vaughan-Preston 
gap is supposed to be located. Far from being evidence 
for the non existence of this feature, this excess must be 
understood as a bias, due to the preferential publication 
of data for objects with certain activity levels, since they 
were gathered from scattered papers of the Mount Wil- 
son group aimed at the study of small samples built from 
different selection criteria. Thus, we cannot join the two 
samples since the representativeness would be lost. This 
problem is very important to our study as long as our ul- 
timate goal is counting stars with certain activity levels, 
after having converted them to ages, to find the star for- 
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Fig. 11. Comparison between the chromospheric activity 
distributions for the S85 and HSDB surveys. 



chromospheric data come from heterogeneous sources, we 
decided to use this sample only as an additional tool for 
our study, in those topics where its inclusion is not likely 
to change the representativeness of the sample. 

In practice, we have to ask whether the initial sample 
itself is representative since to the present moment the 
logi?^ K indices for many stars observed by the Mount 
Wilson group have not been published. We can only base 
our suppositions about the shape of the stellar chromo- 
spheric activity distribution on the two surveys already 
published: S85 and HSDB. The first of these surveys in- 
cludes data for 177 nearby stars in the northern hemi- 
sphere, while HSDB gives data for 817 southern hemi- 
sphere stars. In Figure [ll], we have compared the activity 
distribution of both surveys. Inspection of this plot shows 
that the agreement between these distributions is only fair. 
It is possible to see a separation between active and inac- 
tive stars in both surveys, but the relative number of stars 
in the activity levels seems to be different in them. 

The HSDB survey has selected stars from the combi- 
nation of the two-dimensional MK spectral types in the 



surveys by Houk and collaborators (Houk & Cowley 1975 
~ Houk 119821; Houk & Smith-Moore ~ 



Houk 1978 



with 

the photometric data by Olsen ( 1988| , 1993| ). A secondary 
sample composed by 119 stars, not present in Olsen's pa- 
pers was also observed by them to compensate for the 
incompleteness of these (see below). 

Houk's surveys includes all stars in the Henry Draper 
Catalogue, from which Olsen also has constructed his 
database for photometric observations. The Henry Draper 
Catalogue is supposed to be nearly complete to magni- 
tude V < 9, but Olsen's catalogue has some bia ses to - 
wards more massive stars, as explained in Olsen (1993). 
However, these biases are not expected to depend upon 
any age-related quantity. Thus, the only biases expected 
for our southern hemisphere stars are those present in 
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Fig. 10. Comparison between the main characteristics of the initial (729 stars) and the additional samples (267 stars). 
These main characteristics are the metallicity distribution (upper-left panel), chromospheric activity distribution 
(upper- right panel), and the metallicity-activity distribution (bottom panels). 



ness is V — 8.3 mag, which is the brighter cutoff present 
in Olsen's catalogues (Olsen 1983). For northern hemi- 
sphere stars, the sample is based upon the Catalogue of 
stars within Twenty-five Parsecs of the Sun (Wooley et al. 
197C), which also constitutes a complete sample of nearby 
solar-like stars. 

We are inclined to assume that the HSDB survey best 
represents the real galactic chromospheric activity distri- 
bution, since it includes the largest sample. A comparison 
between the primary and secondary sample in the HSDB 
survey (see their Figure 5), for example, shows that the 
secondary sample (119 stars) present a broader chromo- 
spheric acticity distribution between —4.80 < logi?^ K < 
—5.10, agreeing more with the the S85 survey. However, a 
definitive answer to this question can only be given when 
data for a sample of northern hemisphere stars, as exten- 
sive as that of HSDB, is published. 

Around 87% of our sample is composed of southern 
hemisphere stars, and this makes our chromospheric ac- 
tivity distribution very similar to that of the HSDB sur- 



ours, since we would be comparing very similar samples. 
At least, we can look for some biases if we divide our sam- 
ple and verify whether both subsamples keep the same 
characteristics. We have done this by separating the stars 
according to right ascension: stars with R.A. from 1 ? to 
12 1 ? compose the west sample (410 stars), while that with 
R.A. greater than 12 1 ? compose the east sample (319 stars). 
In Figure [l2] we show the metallicity and log R'n K distri- 
butions and the metallicity-activity diagram for the west 
and east subsamples. Note that the agreement between 
them is much better in all three characteristics we have 
considered. According to this, we can assume that prob- 
ably there is a representative disk chromospheric activity 
distribution. 



2. 7. Spatial Velocities and Orbital Parameters 

For 460 stars from the initial and the additional sample, 
a radial velocity measurement was found in the literature. 
These stars compose the 'kinematical sample'. Each star 
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Fig. 12. Comparison between the main characteristics of the west (R.A < 12 1 ?) and east (R.A > 12 1 ?) subsamples. The 
panels show the same diagrams as in Figure 



the data in the literature, using the equations provided by 
Johnson & Soderblom ( 1987 ). Their orbits were then de- 
termined by numerical integration within a model of the 
Galactic potential, for the calculations of the peri- and 
apo-Galactic distances, R p and R a and the mean Galacto- 
centric radius, R m = (R p +R a )/2 for the orbit (cf. Edv93), 
eccentricity and height above the galactic plane. 

The kincmatical sample is particularly used in a future 
paper of this series aimed to the study of kinematical con- 
straints related to age. Details on this subsample is to be 
given in that papers. 



3. The chromospheric age metallicity relation 

For the derivation of the age-metallicity relation, we have 
to correct the metallicity of the AS and VAS for the mi 
deficiency. Figure |3| shows the age-metallicity diagram 
for the whole initial sample, up to 15 Gyr. The stars seem 
to fall along a very smooth relation just like what is ex- 
pected from chemical evolution theory: the rich stars are 
the youngest, and the poor ones are the oldest. A surpris- 
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Fig. 13. Age-metallicity diagram of the initial sample 
(729 stars). The metallicities of the VAS and the AS are 
corrected for the mj deficiency. Stars with chromospheric 
ages greater than 15 Gyr are not shown in this plot. 

found by Edv93 using isochrone ages. The significance of 
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To find a more refined and unbiased AMR, the sample 
needs further selection. First, we have eliminated 72 stars 
more distant than 80 pc. Taking into account also 11 stars 
that do not have parallaxes measured by HIPPARCOS, 
our sample is reduced to 645 stars. As we wan t to apply a 
volume correction to the AMR (Twarog 1980 ), the sample 
was further limited to apparent V magnitudes lower than 
8.3 mag, by the elimination of 42 stars. 

Eight remaining VAS were also eliminated, since many 
of them may be old close binaries, in which the high ac- 
tivity is produced by highly synchronized rotation, and 
is not related with age. This can explain why some stars 
with ages lower than 0.5 Gyr still present subsolar [Fe/H] 
in Figure [l3]. Rigorously speaking, one may ask why we 
have not discarded the VIS since their very low activ- 
ity might also be unconnected with age and could reflect 
an evolutionary phase analogous to the Maunder Min- 
imum through which the Sun passed during the 17th- 
18th centuries (Baliunas et al. |1995| ; HSDB). We keep 
the VIS in our sample as the effects of the metallicity 
on the log R' HK index make the richer stars resemble older 
IS, or VIS. Many stars amongst the VIS can be normal 
stars, and there is presently no way to separate them from 
the Maunder-mininum stars. The same kind of contami- 
nation by normal stars does not significantly affect the 
VAS grou p, as c an be seen from Figure 6 of Rocha-Pinto 
& Maciel ( |1998|) . 

Seven active stars are known to present greater velocity 
components (Soderblom 1990), and they probably are not 
young stars (Rocha-Pinto, Castilho & Maciel 2000). All 
of these stars were eliminated from the sample. Finally, 
disregarding 37 stars older than 15 Gyr, we arrive at a 
sample with 552 stars. The criteria for elimination are 
summarized in Table |[ 

3.1. Magnitude-limited AMR 

Since our sample is not volume-limited, we have to apply a 
correction to account for the dependence on [Fe/H] of the 
apparent magnitudes. This procedure, called volum e cor- 
rection, was alread y used by some a uthors (Twarog 198C ; 
Ann & Kang 1985 ; Meusinger et al. 1991 ). After binning 
the stars, each metallicity is weighted by d~ 3 , where d is 
the maximum distance at which the star would still have 
an apparent magnitude lower than the magnitude limit, 
which for our sample is 8.3 mag. It is given by 



Table 3. Magnitude-limited AMR. 



8.3-M = 51ogd-5, 



(3) 



where M is the absolute magnitude of the star. 

This magnitude-limited AMR is presented in Table ||, 
where we give the number of stars in each metallicity bin, 
the average- weighted [Fe/H], corrected for the mi defi- 
ciency (see section 2.2), and the metallicity dispersion. 
The metallicity of the last two bins is presented between 
parenthesis, since they are expected to be upper limits 
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Fig. 14. Averaged magnitude-limited AMR (552 stars) 
compared to previous relations in the literature. The 
metallicities of the VAS and the AS are corrected for the 
mi deficiency 



dispersion, we have used the same weights, given by Equa- 
tion (H). These data are presented also in Figure [l4| com- 
pared with previous relations published in the literature. 
In spite of having found a somewhat lower metallicity dis- 
persion, we see a good agreement with the mean points 
of the E dv93's AMR. The agreement with the AMRs by 
Twarog ( 1980 ) and Meusinger et al. ( 1991 ) is marginal. 
These AMRs predict a steady growth of [Fe/H] with time, 
but with a significant flattening in the last 5 Gyr. Instead, 
we have found a steepening in the growth of [Fe/H] at that 
same epoch, a feature also apparent in Edv93. In the old- 
est age bin, our AMR gives an average metallicity of —0.48 
dex which is 0.14 dex higher than the average metallicity 
in Edv93. We believe that this discrepancy results from 
the errors in the chromospheric ages (see discussion on 
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Criterion stars eliminated 



Stars not having parallaxes in the HIPPARCOS database 11 

Stars distant by more than 80 pc 72 

Apparent V magnitude greater than 8.3 mag 42 

Very active stars (log R' HK > -4.20) 8 

Objects known to be chromospherically young, but kinematically old 7 

Stars with nominal chromospheric ages greater than 15 Gyr 37 



We have not found an absence of metal-rich stars with 



ages between 3 and 5 Gyr (Carraro et al. 1998 ). In fact, our 
metal-rich stars are very concentrated at the younger bins, 
since this AMR does not show the large scatter present in 
other studies. 

3. 2. The initial metallicity of the disk 

The AMR presented in Table || indicates a high estimate 
for the initial metallicity of the disk. According to it, we 
should expect an average [Fe/H] rj —0.48 dex at around 
13.5 Gyr ago. This metallicity is somewhat higher than the 
corresponding values found by Twarog (1980) and Edv93. 

A high initial metallicity is indicative of significant pre- 
enrichment of the gas before the formation of the first stars 
in the disk. That the disk has had some previous enrich- 
ment can be easily seen from the lower cuttoff in its metal- 
licity distribution in [Fc/H] ss —0.8 dex (Rocha-Pinto & 
Maciel 1996). Even in the framework of an infall model, 
the disk initial metallicity must be non-zero in order to 
match the G-dwarf problem. 

Our AMR has a peculiar behaviour towards greater 
ages. It flattens while the other relations (as well as the 
theoretical models) generally becomes steeper, indicat- 
ing a rapid enrichment in the early galactic phases. That 
flatenning is the major difference between the mean points 
of our AMR and the mean points of that by Edv93. 

The behaviour of our AMR in the oldest age bins are 
probably very affected by the age errors discussed in sec- 
tion 2. The chromospheric activity yields more accurate 
ages for younger stars, and we expect that the AMR at 
the youngest bins is fairly well reproduced. For the older 
bins, however, we have to estimate the mean deviation 
from the real AMR that we expect to find from using 
chromospheric ages. This can be done by simulating the 
scattering of the stars in the AMR due to the age errors. 

The procedure follows closely that used in Paper II to 
compute the statistical confidence levels for the star for- 
mation history (hereafter SFH) features, and we refer to 
that paper for more detail. We simulate a constant SFH 
composed by 3000 stars. For each star, a metallicity is 
assigned by a pre-adopted AMR. After that, the stellar 
age and its metallicity are used to derive the correspond- 




Age (Gyr) 

Fig. 15. Simulation of the AMR, with a set of 700 stars. 
The solid curve shows the original linear AMR. The points 
with error bars indicate the recovered AMR. It can be seen 
that the chromospheric age errors introduce a small scat- 
ter in the AMR, but do not destroy the average relation. 

Maunder-minimum phase. This is done by inverting the 
equations presented by Rocha-Pinto & Maciel ( 1998| ). A 
database composed by 3000 stars with logR' HK and [Fe/H] 
is then built. The stars in this database are binned in 0.2 
Gyr intervals, and we discard randomly in each interval 
the number of stars expected to have died or to have left 
the galactic plane at that corresponding age (see Paper II). 
The remaining stars compose a sample of around 600-700 
stars. For these stars, we randomly shift their correspond- 
ing logR' HK and [Fe/H] according to the corresponding 
errors in these quantities. The final catalogue then resem- 
bles very much the real data sample used in this study. 

The final simulated catalogue is used to derive the 
AMR in the same form we did for the real data. For the 
sake of simplicity, we assume a simple linear law for the 
AMR: 



[Fe/H] = 0.44-0. lOt. 



(4) 



where t is the stellar age in Gyr. 

The simulations were repeated 20 times. The results 
are very similar from one simulation to the other, and one 
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that the AMR found with the chromospheric ages follows 
closely the real AMR of the disk up to 9 Gyr ago, when 
it begins to deviate. The deviation is always in the sense 
of an increased dispersion and higher mean AMR. This is 
exactly the same that is observed in the AMR we have 
found for the disk. On the other hand, these simulations 
show that even in face of great age errors, the AMR can 
be fairly well recovered from a sample with chromospheric 
ages. 

The increase of metallicity dispersion at the older age 
bins, as well as the greater deviation from the real AMR, 
are all consequences of the age errors. The older bins tend 
to be populated by the poorer stars, and the corrections 



to the chromospheric ages (Rocha-Pinto & Maciel 1998) 
increase very much for metal-poor stars. A small error in 
[Fe/H] or logi?^ K reflects in a substantial error in age, 
which most probably will push the star to greater chro- 
mospheric ages. Due to this effect, the oldest age bins are 
more likely to be depopulated by metal-poor stars than 
the other bins. That is why the recovered AMR fails to 
match the original AMR as we consider ages progressively 
older. And the metallicity dispersion increase because the 
absolute age errors increases as a function of age. 

In our simulations we have found that the average 
metallicity of the oldest bin is invariably around 0.20-0.23 
dex higher than the real average metallicity, due to the 
age errors. If we apply this value to the Milky Way AMR 
derived previously, we get a disk initial metallicity around 
—0.70 dex, which agrees very well with Edv93's AMR and 
with the G dwarf metallicity distribution. 



4. Metallicity dispersion 

4-1- Some Evidences for a Non- Homogeneous Interstellar 
Medium 

The metallicity dispersion we have found is about 0.13 
dex as shown by the error bars of Figure H and fr om 
Table [| It is much lower than that found by Edv93. At 
first sight, this result seems to revitalize old ideas about 
the chemical evolution of our Galaxy, in which the inter- 
stellar medium has been continuously and homogeneously 
enriched by metals ejected by the stars. Previous works 



on the AMR (Twarog [198$ Meusinger et al. |l99l| ) have 
consolidated this view by finding crpo/H] around 0.12-0.18 
dex. 

This picture was strongly questioned after Edv93 pub- 
lished their detailed work on the chemical evolution of the 
disk. They found a larger dispersion in the AMR, varying 
from 0.18 to 0.26 dex. A number of hypotheses to explain 
it were considered, and the main conclusion of the authors 
is that a significant part of this scatter could be physical. 
Infall is quoted as one of the best mechanisms that could 
drive such an intrinsically large metallicity dispersion at 
all epochs, if the time scale for the mixing of the infalling 
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spectroscopic [Fe/H] 

Fig. 16. Comparison between photometric and spectro- 
scopic metallicities for Edv93's stars. The photometric 
metallicity for them was estimated as described in sub- 
section 2.2. 



the infall hypothesis requires star formation to be much 
more efficient in the material that has just fallen onto the 
disk than in the already well- mixed gas. This would oc- 
cur if infall could drive star formation. It is interesting 
to see that a number of theoretical and observational evi- 



dences favour this process (Lepine & Duvert 1994 ; Lepine 
et al. 1999| ), although a detailed study of this star for- 
mation mechanism connected with a chemical evolution 
model has never been made. 

An alternative explanation was proposed by Wielen et 
al. ( 1996 ) and Wielen & Wilson ( 1997 ) according to which 
the metallicity dispersion in Edv93's AMR origin ated from 
the diffusion of the stellar orbits (Wielen 1977 ). Accord- 
ing to this hypothesis, the Sun was born 1.9 kpc closer to 
the galactic center in comparison with its present position. 
This could solve a long-lasting puzzle about the fact that 
the Sun is richer than some younger neighbour objects 
(Cu nha fc Lambert |1992| ; de Freitas Pacheco [l993| ) . Clay- 
ton (1997) also used this scenario to explain the meteoritic 
abundance ratios of the isotopes 29 Si and 30 Si, related to 
28 Si. 



Recently, Binney & Sellwood ( |2000[ ) reconsidered the 
diffusion of stellar orbits and found that, a typical star is 
unlikely to migrate from the galactocentric radius of its 
birthplace by more than 5% over its lifetime. These au- 
thors say that Wielen et al. ( 1996| ) calculated the diffusion 
of stellar orbits in velocity space, adopting isotropic and 
constant diffusion coefficients. When the diffusion is cal- 
culated in integral space, it becomes very anisotropic. As 
a result, the star does not change significantly the guiding- 
center radius of its orbit, which should be similar to R m , 
as in Edv93. This recent work concludes that the metal- 
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The existence of some scatter in the interstellar 
medium is not questioned. The Sun-Orion abundance dis- 
crepancy and the existence of young stars with subsolar 
abundance (Grigsby, Mulliss & Baer 1996 ) are classical 
puzzles that point to the existence of an intrinsic metal- 
licity dispersion in the galactic gas. The real problem is 
its quantification, since these anomalies can be outliers. 

For instance, Binney & Sellwood ( 2000 ) show data 
measured by other authors, which indicate that the intrin- 
sic scatter in [ O/H] should be around 0.10 dex. Garnett 
& Kobulnicky ( [2000 ) also conclude that, from measure- 
ments in near by galaxies and the local ISM (Kennicutt & 
Garnett 1996; K obulnicky & Skillman 1996; Meyer, Jura 



& Cardelli 1998 ), the metallicity scatter is lower than 0.15 
dex. 



4-2. Isochrone- spectroscopic data vs. 
chromo spheric- photometric data 

The average metallicity dispersion we ha ve found is very 
close to that found by Twarog (1980; ~ 0.12 dex). 



Twarog's sample is greater than Edv93's and has a good 
statistical significance, but his metallicities are less accu- 
rate. Moreover, his ages were found by old, now outdated, 
isochrones. The problem can be summarized as follows: 
photometric AMRs show a small, well-behaved metallic- 
ity scatter (however, see Marsakov et al. 1990), while the 
only one based on a spectroscopic sample indicate the op- 
posite. Note that Ng & Bertelli ( p.998| ) only revisit Edv93's 
ages, so that their AMR cannot be taken as an indepen- 
dent evidence for a real greater scatter. 

It is important to show that there is no significant 
difference in the quality of our data compared to those 
of Edv93, even taking into account that our metallicities 
are photometric and our ages are chromospheric. In Fig- 
ure |l^, we show the comparison between the 16 stars in 
common between ours and Edv93's sample, before and 
after the application of the metallicity corrections to the 
chromospheric ages (Rocha-Pinto & Maciel 1998). Note 
that these corrections improve substantially the matching 
between both methods to measure the stellar ages. The 
exceptions are few, and will be discussed separately. 

The agreement is also good for metallicities. In Figure 
|l8| , the position of our stars in the age-metallicity diagram 
is shown. Lines connect these stars with their position in 
the same diagram using Edv93's data. With only one ex- 
ception, all stars present differences in metallicities that 
are within the expected error in the photometric calibra- 
tion of [Fc/H]. 

Sixteen stars is a small number to test if our method 
to estimate ages and [Fe/H] is good. We have chosen to 
compare directly these methods to estimate stellar ages 
and metallicities, using the stars in Edv93's sample. 

The first of these comparisons is shown in Figure [h| 
where we show the correlation of photometric and spectro- 
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Fig. 18. Age-metallicity diagram showing the position of 
the 16 stars in common between this work and Edv93. The 
lines connect the same stars as they appear using photo- 
metric plus chromospheric data and isochrone plus spec- 
troscopic data. With only one exception, all stars present 
differences in metallicities that are within the expected 
error in the photometric calibration of [Fe/H]. 



photometric metallicity was estimated as described in the 
subsection 2.2. The agreement is typical of that expected 
from a photometric calibration. The standard deviation of 
the data is 0.10 dex. 

Chromospheric indices were published only for 81 stars 
from those of Evd93 sample, and 40 of them were used in 
the determination of the metallicity cor rectio n to the chro- 
mospheric age (Rocha-Pinto & Maciel 1998 ). We present, 
in Figure [l^, a comparison between the isochrone and the 
chromospheric ages, before and after the application of 
the metallicity corrections. The stars are distinguished 
by symbol (according to their spectroscopic metallicity) 
and style (if it was already used by Rocha-Pinto & Maciel 



1998 or not). The Figure shows clearly that the metallicity 



corrections improve substantially the chromospheric ages. 
The scatter is comparable to the expected error in both 
methods (we use here a formal average error of 0.1 dex for 
Edv93 ages, although Lachaume et al. 1999 have shown 
that this error can be easily underestimated). Panel c of 
Figure |l9| sho ws the same comparison, now using Ng & 
Bertelli (1998) ages. The agreement is somewhat better, 
although a greater isochrone age is still found for some 
stars. Panel d shows a comparison, in the same scale, be- 
tween the isochrone age determinations by Edv93 and Ng 
& Bertelli ( gjjg§ ). 

Isochrone ages are better compared to chromospheric 



ages, at least for early G dwarfs (Lachaume et al. 1999 



showed that this is not true for intermediate and late G 
darfs). The difference can be large, but in general, it is of 
the same order as the difference between two isochrone age 
determination made using different isochrones (see panel 
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Fig. 17. Comparison between isochrone and chromosphcric ages for 16 stars in common between this work a nd Ed v93. 
a) before the application of the metallicity corrections to the chromosphcric ages (Rocha-Pinto & Maciel 1998 ); b) 
after the application of the metallicity corrections. The few deviating stars are discussed in section 4.3. 



in Fig. H is fairly reasonable. They are considerably large, 
but they are not expected to destroy the AMR, as shown 
by our simulation in Figure |lj. 

There are only three stars that were classified as sub- 
giants in Figure p], that have both a chromospheric and 
isochrone age. For two of them (HD 9562 and HD 131117), 
the ages agree remarkably to within 0.15 dex. The third 
star (HR 4734 = HD 108309) is investigated in more de- 
tail in what follows. On the other hand, Edv93's sample 
has some slighlty evolved subgiants, which are represented 
by the oldest stars in Figure The good agreement be- 
tween their chromospheric and isochrone ages shows that 
the chromosphcric activity-age relation can be extrapo- 
lated to slightly evolved subgiants. 

4-3. Anomalies to the chromospheric activity-age relation 

Few stars deviated significantly from the expected rela- 
tion in panels b and c of Figure |l^. With two exceptions 
(HR 1780 and HD 165401), all of them have metallici- 
ties greater than —0.05 dex. The deviation is in the same 
direction: stars that appear to be very young, according 
to their chromospheres, are found to be much older from 
their position in the HR diagram. Three of these stars (HR 



by the Mount Wilson group, according to Duncan et al. 
( 1991 ). While their high chomospheric activity could be 
caused by a periodic active phase in their magnetic cycles, 
all of the other stars were observed more than a hundred 
times, over a time span of 12 years or more, and errors 
in the chromosphcric indices cannot be invoked to explain 
their different ages. 

Some of the deviating stars in panel b are chromo- 
spherically young, kinematically old stars (hereafter re- 
ferred as CYKOS), that have been recently investigated by 
some of us (Rocha-Pinto, Castilho & Mac iel pOOOp . Their 
origin is presently unknown. Poveda et al. ( 1996| ) have also 
found objects like these amongst UV Ceti stars, and they 
proposed they could be low-mass blue stragglers, formed 
by the coalescence of close binaries. Other explanation 
could be that these stars are themselves close, unresolved 
binaries. The high rotation rates of these systems tend to 
keep the chromospheric activity well beyond the normal 
age where we would expect it to have been diminished. 
Moreover, systematic chemical variations between them 
and the other normal stars could be a cause for some de- 
viations from the relation. However, we have found no such 
difference, using the abundance ratios provided by Edv93. 
We can but conclude that these stars do not follow the 





log Isochrone Age log Age (Edvardsson et al.) 



Fig. 19. Comparison between chromospheric and isocrone ages for 81 Edv93's stars. The different symbols stand for: 
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explored by future investigations is whether the metallic- 
ity correction to the chromospheric ages of the metal-rich 
stars have been overestimated. It is apparent from Figure 1 
of Rocha-Pinto & Maciel ( |1998|) that some metal-rich stars 
present great difference between their chromospheric and 
isochrone ages. 

The existence of anomalies to the chromospheric 
activity-age calibration does not throw doubt about the 
use of chromospheric indices to measure stellar ages. The 
exceptions are few. There are also exceptions to the use 
of photometric indices to measure stellar parameters, for 
instance in peculiar stars, but these do not rule out the 
validity and quality of photometrically calculated param- 
eters. The difference is that there are ways to photomet- 
rically distinguish between a normal from a peculiar star. 
But this is not possible for CYKOS from chromospheric 
indices alone. We are not sure whether it could be done by 
photometry. The only way to discover a CYKOS presently 
is by comparing their activity levels with their isochrone 
ages and their kinematical characteristics. 

4-4- The problem of the metallicity scatter 

A more instructive way to compare the metallicity scat- 
ter of the two AMR is presented in Figure |2(], where we 
superpose both AMR's. The differences occur mainly in 
two regions, which we have marked by dotted lines, and 
named as Region I and II. Region I has also other interest- 
ing characteristic: it is well detached from the bulk of stars 
in both samples. Note that four stars amongst those that 
have deviated from the mean relation in Figure [l^ are also 
present in Region I. It suggests that the other stars in this 
location do not follow the chromospheric activity-age re- 
lation. We have chromospheric indices for five other stars 
in this region (HR 3176, HR 3951, HR 5423, HR 8041 and 
HR 8729). All of them also present chromospheric ages 
lower than their isochrone ages, although the age excess 
is smaller than that found for the stars in Figure [l9]. Ta- 
ble U presents the logarithmic isochrone age of these stars 
(log/age), in Gyr, taken from Edv93, and the age excess 
(A), which we define as the logarithmic difference between 
the isochrone and the chromospheric age of a star (see 



Table 4. Age excesses for stars in Region I and II of 
Edv93's AMR. 



Rocha-Pinto & Maciel |1998| ). For other star, HR 2354, 
chromospheric indices were not found in the literature, 
but emissi on in the Ca II K line was observed by Hagen 
& Stencel ( 1985 ). Significant X-ray emission was also ob- 
served bythe ROSAT observatory (Hiinsch, Schmitt & 
Voges 1998j ), confirming that it is chromospherically ac- 
tive. Its spectr al typ e is G3 III-IV. Although Micela, Mag- 
gio & Vaiana ( 1992] ) showed that the X-ray activity can 
be used as an age indicator in giants, no study has ever 
aimed to calibrate an age relation for them. We have u sed 
the relation proposed by Kunte, Rao & Vahia ( 1988 ) to 
estimate an X-ray age for it. We alert that this relation 
was calibrated only for main-sequence stars, and its use 



Name 


log J agc 


A 




Ftcgion I 




HR 2354° 


0.84 


0.44 


HR 3176 


0.94 


0.31 


HR 3538 


0.95 


0.87 


HR 3951 


0.81 


0.48 


HR 4734 


0.93 


0.59 


HR 5019 


1.15 


0.55 


HR 5423 6 


0.84 


0.53 


HR 7232 


1.08 


0.79 


HR 8041 6 


0.89 


0.48 


HR 8729 


0.93 


0.53 




Region II 




HR 4657 


0.73 


0.32 


HR 5447 


0.32 


-0.17 


HR 8354 


0.75 


-0.32 



a Age estimated by X-ray luminosity (Hiinsch e t al 
b Chromospheric indices given by Young et al. ( 1989 ) 



1998) 



luminosity over bolometric luminosity given by (Hiinsch, 
Schmitt & Voges 1998 ), we have an X-ray age of 2.5 Gyr 
to HR 2354, which is lower than the isochrone age found 
by Edv93, just like for the other stars with chromospheric 
ages. 

The presence of CYKOS in this region also points to 
a wrong age determination (both isochrone and chromo- 
spheric), since as a coalesced star, as an unresolved close 
binary would have both chromospheric activity and po- 
sition in the HR diagram uncorrelated with its real age. 
This can explain why there is a gap between Region I and 
the other stars. Chromospheric indices for other stars in 
this region could test more properly this hypothesis. 

The other region marked in Figure ^0| corresponds to 
the metal-poor stars with young to intermediate ages. 
None of such stars are present in our AMR. This could be 
explained if there were some systematic error in the metal- 
licity corrections to the chromospheric ages that would 
avoid locating stars in this region. Unfortunately, there 
are chromospheric measurements for only three stars in 
Region II. Their age excesses do not follow a systematic 
trend, as for the stars in Region I, as can be verified in 
Table |. 

We have checked if the scatter in our AMR could have 
been partially destroyed due to the use o f the metallicity 
corrections. Rocha-Pinto & Maciel ( 199£ ) have proposed 
such a correction from their finding that the difference 
between chromospheric and isochrone ages were related 
to the stellar metallicity. The correction itself disregards 
the scatter in the data. If this scatter is real, probably re- 
flecting different chemical compositions, the use of a cal- 
ibrated metallicity correction could remove partially the 
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simulation analogous to that used in the previous section. 
This time, we have used a 'AMR' composed only by scat- 
ter, that is, 730 stars were randomly distributed in age 
(from to 16 Gyr) and metallicity (from —1.2 to 0.4 dex). 
We have taken explicitly into account the dispersion in 
the metallicity correc tion, taken from Figure 1 of Rocha- 
Pinto & Maciel ( 1998| ), as well as the errors in log i?HK an d 
[Fe/H]. We have found that the chromospheric AMR pre- 
serves closely the real metallicity dispersion of the data, 
even with the greatest possible scatter for the disk AMR. 

We have found no indication that the differences in the 
AMR are due to difference in the methods used to estimate 
age and metallicity. We are left with the hypothesis that 
one of the samples (Edv93's or ours) is not suitable to 
investigate this constraint. 

There are arguments to suspect that Edv93's sample 
selection has biased their results. Their sample was se- 
lected very carefully to allow the study of the abundance 
ratios evolution. It is not suitable to other studies such as, 
for instance, the metallicity distribution. Part of the scat- 
ter in Edv93's AMR comes from their selection procedure 
to have nearly equal number of stars in predetermined 
Smi intervals, to assure a good coverage of [Fe/H] values 
in their sample. The authors themselves make a correction 
for this, using a metallicity distribution of 446 dwarfs, af- 
ter which the metallicity dispersion diminishes to around 
0.21 dex. While this value is closer to previous estimates 
of the metallicity dispersion in the solar vicinity, it is still 
greater than ours. What is generally not recognized is that 
the same kind of problem can occur in the investigation 
of the metallicity dispersion inside each age bin. 

We have tested this with a number of artificial samples. 
We have deliberately selected randomly nearly 25 stars in 
each of the metallicity bins that Edv93 used, amounting 
189 stars just like their sample, from a set of 2000 artificial 
stars. These pre-selected AMRs were compared with other 
AMRs composed by 189 stars randomly selected from the 
same parent population. We have verified that contrary to 
the expectation, this kind of selection procedure does not 
necessarily imply a metallicity scatter greater than what 
would be found with a more unbiased selection procedure. 
However, there is a significant difference between these 
simulations and the procedure made by Edv93. Accord- 
ing to them, they had observed the 25 brightest stars in 
each metallicity box, after disregarding binaries and stars 
rotating faster than 25 km/s. This probably causes a bias 
that was not taken into account by the authors. If there 
was a one-to-one relation between metallicity and age, the 
age of the 25 brightest stars would be the same as that of 
the 25 faintest stars. On the other hand, if there is a real 
metallicity scatter in the interstellar medium, at the same 
metallicity, the brightest stars would not have the same 
age dispersion as the faintest stars, or the age dispersion of 
a subsample randomly selected. The brightest stars will be 
generally the nearer stars and the most evolved, which are 



ory we also know that these are generally the older stars, 
or the hotter stars. Since Edv93 also avoid faster rotating 
stars (which corresponds to the hotter stars), the majority 
of these stars will be old. Therefore, Edv93's metallicity 
bins have probably non-representative age distributions, 
which changes the whole dispersion in the AMR. 

This selection effect was introduced by intrinsic ob- 
servational limitations. Other selection effects are also 
present in their sample due to this problem. For instance, 
Garnett & Kobulnicky ( 2000| ) have shown that there is a 
distance bias in Edv93's sample. Stars at distances 30-80 
pc from the Sun are systematically metal-poor. The expla- 
nation for this trend is provided by Edv93 themselves. To 
have equal number of stars in each metallicity bin, they 
observed the metal-poor stars at a large volume. Garnett 
& Kobulnick conclude that the real scatter in the AMR 
should be lower than 0.15 dex, which favours our results. 



Edmunds 1995 favours the same conclusion. 



However, the puzzle cannot be considered solved. It is 
highly necessary to have an independent study of the spec- 
troscopic AMR, preferably using a sample as large as that 
of Twarog. It is also important to make a more detailed 
investigation of the link between age and stellar activ- 
ity, especially for the CYKOS and other deviating stars. 
Presently, we can conclude that the photometric, chromo- 
spheric AMR does not give support to the great scatter in 
the AMR as found by Edv93. Therefore, the AMR can be 
still regarded as a strong constraint to chemical evolution 
models, just like the G dwarf metallicity distribution and 
the radial abundance gradients (see Chiappini, Matteucci 
& Gratton|l997|). 



5. Conclusions 

We have used an extended sample composed of 552 late- 
type dwarfs, with chromospheric ages and photometric 
metallicities, to address the AMR of the solar neighbour- 
hood. Our main conclusions can be summarized as follows: 

1. The AMR is found to be a smooth function of time. 
The average metallicity has increased by at least 0.56 
dex in the last 12-15 Gyr. No absence of rich stars age- 
ing 3 to 5 Gyr, as suggested by Carraro et al. (1998), 
was found. The results are in fairly good agreement 
with t he mean points of Edv93's and Ng & Bertelli 
(tL998|) 's AMR. 

2. The initial metallicity of the disk was around —0.70 
dex. This suggests some previous enrichment of the 
gas which gave rise to the thin disk. 

3. A number of very metal rich stars in Edv93 could 
be composed by chromospherically young, kinemati- 
cally old stars (CYKOS). These stars could originate 
from coalesced stars or unresolved close binaries. This 
may indicate that their isochrone ages, as measured 
by Edv93, are wrong. This can, at least, explain why 
Edv93 have found some old stars more metal rich than 
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Fig. 20. AMRs from Edv93 (crosses) and from this work (open circles). The major differences correspond to two areas 
in this diagram where we have found no stars, while Edv93 found some. These regions were named I and II, and are 
marked by dotted lines. The stars labelled have also shown peculiar behaviour in Figure 111]. 



4. An average intrinsic metallicity dispersion of 0.13 dex 
was found. Several hypotheses were tested to explain 
this scatter compared to that found by Edv93. We 
have found no indication that it could be produced by 
our method to estimate ages and metallicities. On the 
other hand, this dispersion agrees closely with the re- 
cent conclusions by Garnett & Kobulnicky ( 2000| ) that 
the real scatter in the AMR should be lower than 0.15 
dex. Nevertheless, additional independent determina- 
tions of the AMR are strongly encouraged to confirm 
this view. 

Some important bias in the sample can be investigated 
with the aid of kinematical data. A third paper on this 
series addressing this topic is being planned. 
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